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We argue that the magnetic susceptibility data, Refs.— for the low-density two-dimensional (2D) 
silicon-based electron gas indicate that magnetically active electrons are localised in spin-droplets. 
The droplets exist in both the insulating and metallic phases, and interact ferromagnetically, forming 
an effective 2D Heisenberg ferromagnet. Comparing the data with known analytical and numerical 
results for a 2D Heisenberg ferromagnet, we determine that JS 2 fa 0.6 K, where S is the spin of the 
droplet and J is the ferromagnetic exchange constant between droplets. We further argue that most 
likely S = 1 with four electrons occupying each droplet on average. We discuss the dependence of 
the magnetic susceptibility and the specific heat on the external magnetic field, which follows from 
the model, and hence we suggest further experimental tests of the model. 
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I. INTRODUCTION 

Studies of low density two-dimensional electron gas 
(2DEG) systems attract great attention, because of the 
unusual and rich properties of these systems. In the 
present work, we consider the magnetic properties of a 
silicon based 2DEG discovered in recent studies 1 - - — . Most 
likely the properties are related to the metal-insulator 
phase transition (MIT), and, in our opinion, understand- 
ing of these properties sheds some light on the nature of 
the transition. 

The MIT in a 2DEG has attracted protracted atten- 
tion from both experiment and theory, and remains a 
puzzling area of research to date^~— . It was once believed 
that a MIT in such systems could not take place, because 
a true metallic phase does not exist in a non-interacting 
2DEG § , although extension of the scaling theory of local- 
isation to include the effects of interaction^^ suggested 
that a MIT may be possible. Resistivity measurements 
in a silicon-based 2DEG finally provided evidence for a 
true MIT—, with numerous works following thereafter 
and studying the transition-^—. The mechanism for a 
MIT in a 2DEG remains unclear to this date, but the ex- 
istence of localised states on the insulator side of a MIT 
is generally accepted£&2£~— . 

Intimately linked to the problem of the MIT is the na- 
ture of the ground state of a 2DEG, which still remains an 
outstanding proble m 20 : 22 " —. The ground state depends 
on electron density. There is little doubt that at a suf- 
ficiently high density it is a normal paramagnetic Fermi 
liquid 2 ^. At a lower density, the system might have a 
Stoner transition to a ferromagnetic Fermi liquid 2 ^, and 
ultimately at a very low density it must undergo a tran- 
sition to the Wigner crystal^ .. These are scenarios for 
a 2DEG without any extrinsic disorder, see Refs<2i~— . 
Extrinsic disorder can further complicate the situation, 
see, for instance, Rets . 35 ' 36 . 

The magnetisation and magnetic susceptibility of a sil- 
icon based 2DEG in the electron density range below and 
above a MIT have been studied recently^— The exper- 



iments have been performed with an in-plane magnetic 
field and so only spin related magnetic properties have 
been measured. There are three important outcomes of 
these measurements, (i) Thermodynamic magnetic prop- 
erties vary continuously across the MIT. (ii) The zero- 
field magnetic susceptibility diverges rapidly in the limit 
T — > 0, xo oc 1/T 2 - 4 , in both the insulating and metallic 
phases, (hi) In the metallic phase at T ~ 1 K, the value 
of the susceptibility is by orders of magnitude larger than 
the expected value of the ideal gas Pauli susceptibility. 

The authors of Refsii - — explain their data through 
the formation of electron droplets. Each droplet has a 
nonzero spin. These droplets melt in the metallic phase 
with increasing density and temperature, but continue to 
exist up to large densities. At a fixed density of droplets, 
this picture would give the usual Curie scaling of the sus- 
ceptibility with temperature, oc 1/T. To explain the ob- 
served scaling, oc 1/T 2 4 , Ref^ suggests that the density 
of droplets is decreasing when temperature is increasing. 

In the present work, we take a somewhat different view 
to explain the magnetic data. We agree that the data 
practically unambiguously indicate the formation of elec- 
tron droplets with nonzero spin. The droplets are proba- 
bly formed due to extrinsic disorder or extrinsic disorder 
assisted by the Coulomb interaction. Pragmatically for 
our purposes, the exact mechanism of their formation is 
not important. We only assume that at low temperatures 
all the internal degrees of freedom of the droplets are 
frozen and hence the only dynamical degree of freedom 
is the spin of the droplet. The very steep temperature 
dependence of the magnetic susceptibility, in our opin- 
ion, indicates ferromagnetic instability. So we assume 
a ferromagnetic interaction between droplets and hence 
consider the system as a quantum Heisenberg model on a 
random 2D lattice. This is a very natural model for the 
insulating phase and we also assume that the localised 
droplets described by the Heisenberg model continue to 
exist in the metallic phase. While in the metallic phase 
most electrons go into itinerant states, these electrons are 
magnetically almost idle and the main contribution to 
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the magnetic susceptibility is due to the relatively small 
fraction of electrons localised in the droplets. The den- 
sity of the droplets diminishes in the metallic phase with 
increasing electron density. 

The structure of the paper is as follows. In Section II, 
we review known properties of the 2D quantum Heisen- 
berg ferromagnet and derive some previously unknown 
properties, which we need for the present work. A special 
point of the 2D Hcisenbcrg ferromagnet is exponential de- 
pendence of the susceptibility on a special combination of 
spin and effective exchange constant. Comparison with 
experimental data performed in Section III allows us to 
determine this combination rather accurately. The num- 
ber of electrons per droplet and the spin of the droplet 
appear only in a prefactor to the susceptibility. There- 
fore, we provide additional arguments to determine their 
values. In Section IV, we make predictions, which can be 
checked experimentally. Finally, Section V presents our 
conclusions. 



II. 2D QUANTUM HEISENBERG 
FERROMAGNET 

The Heisenberg model is defined by the Hamiltonian 
Hj = - J Si ■ Sj . (1) 

<ij> 

We assume that the model is defined on a square lattice. 
Summation in (TJJ) is performed over nearest sites and Si 
is the quantum spin at the site i. In the ground state, all 
spins are aligned ferromagnetically along, say, the z axis. 
Excitations are spin waves with the following spectrum^ 

e k = 2JS [2 - cos(/c x ) - cos(fc y )] JSk 2 . (2) 

Hereafter we set the Planck constant and the Boltzmann 
constant equal to unity, K = ks = 1. 

Each spinwave excitation carries spin AS Z = — 1. The 
excitations are bosons and therefore the magnetisation 
at a nonzero temperature is 

{Sz)=S ~ I e^-l(2ny (3) 

The integral is logarithmically diverging at small mo- 
menta. This is a direct consequence of the Mermin- 
Wagner theorem^, which claims that long range order 
is impossible in a 2D system at a nonzero tempera- 
ture. The ferromagnetic ordering exists only within a 
correlation length £. To find value of £, one has to set 
(S z ) = and impose a lower limit in the integration in 
k > k m i n ~ This gives the following correlation 
length in the low temperature limits 

i oc e 2 * Js2 ' T . (4) 



There are N ~ £ 2 spins within the correlation length, 
these spins act as a magnetic domain with total mag- 
netic moment M ~ SN. The concentration of do- 
mains is no ~ 1/N. All in all, this describes a super- 
paramagnet with the following magnetic susceptibility, 
X oc ri£,M 2 oc N oc £ 2 oc e 47T,JS l T . This simple logic 
does not give the prefactor before the exponential. The 
renormalisation group (RG) calculation gives^ 

where A is a constant. Note that the third power of 
the semiclassical parameter 47r j g a in the prefactor in 
([5]) arises in the two-loop approximation, while the sin- 
gle loop approximation gives only the first power of the 
parameter—. 

Quantum Monte Carlo (QMC) simulation of the mag- 
netic susceptibility for 5 = 1/2 was performed in Ref.— 
and for S = 1 in Ref4i In our analysis, we use the results 
of Rcf4i, because this simulation accounts for a nonzero 
magnetic field. The susceptibility QMC data^L for the 
values of magnetic field B = 0.005J and B = 0.01J are 
presented in Fig. [TJ 




FIG. 1: Magnetic susceptibility per site obtained in QMC 
simulations of S — 1 Heisenberg model^. Values of the mag- 
netic field are B — 0.005J (red circles) and B = 0.01J (black 
squares). The green dashed line shows fit to the RG Eq. ([5|. 
The solid blue line shows the simple exponential fit (|10[) . 

Interaction with a magnetic field is defined by the 
Hamiltonian 

H B = Hj-Bj2Siz ■ (6) 

i 

Even a very small magnetic field significantly influ- 
ences the susceptibility at low temperature. Eq. ([5]) is 
valid only at low temperature, JS > T, but on the other 
hand, due to the presence of a magnetic field, the tem- 
perature cannot be too low, T > T m . Here T m is the 
temperature where the susceptibility is maximum, see 
Fig. [TJ The value of T m depends on the magnetic field. 
The very strong dependence of the susceptibility on the 
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magnetic field is related to the dimensionality of the sys- 
tem. To explain the dependence, we remind that at zero 
temperature the magnon dispersion in a magnetic field 
at small k is2I 



e k = JSk 2 



B. 



(7) 



When deriving Eqs. (j4} and (O, we substitute in ([3]) 
kmin ~ l/£ as the lower limit of integration. This is 
correct only if JS/£ 2 > B. From this condition, one 
immediately finds T m 



AttJS 



hi 



' i-n.JS 2 ' 
* B , 



(8) 



The dependence on the magnetic field is logarithmic and 
hence even a tiny magnetic field significantly influences 
the susceptibility. The authors of Ref4i used a power 
function to fit the dependence of T m on the magnetic 
field, T m oc B 1 . The fit gave a pretty small value 
7 ps 0.15. While at B ~ J, the power fit might make 
sense, at small fields, B <C J, the dependence is certainly 
logarithmic. 

The RG expression (|5|) for magnetic susceptibility is 
valid at low temperature, T -C JS, as it follows from the 
dispersion ((2|) . On the other hand, the temperature must 
be higher than T m . So the region of validity of Eq. ([5} is 



T m < T < JS 



(9) 



Clearly, the right slope (green dashed line) in Fig. [T] is 
not quite in this region, the magnetic field in these QMC 
data is not sufficiently small and hence the temperature 
is not sufficiently low. Nevertheless, we try to fit the data 
using Eq. ([5]). The RG fit shown in Fig. [T]by the green 
dashed line is not bad. However, the simple exponential 
fit 
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0M2Se 7 - 6Js2/T 



(10) 



shown by the blue solid line is better. Below in the anal- 
ysis of experimental data we will use the simple exponen- 
tial fit (|10p . having in mind that the data are also taken 
at not very low temperatures. Practically, this means 
that we will compare the experimental data directly with 
results of QMC simulations. Note that Eq. (fTU|) at 
T 3> JS 2 is approaching a constant instead of Curie's 
law, x = S(S + 1)/(3T). Hence Eq. (JTDJ) overestimates 
the susceptibility at large T. In Eq. (fTU|) , we assume lin- 
ear scaling of the prefactor with S as it is predicted by 
RG approach, see Eq. ([5]). This is the scaling at a fixed 
value of the semiclassical parameter JS 2 /T. 

The electron spin droplets in silicon certainly do not 
form a regular square lattice. However, in the temper- 
ature range 4nJS 2 3> T > T m , the exact structure of 
the lattice is not important. The dynamics depend only 
on the quadratic magnon dispersion, = JSk 2 , where 
J is an effective exchange constant. The dispersion is 
the universal property of a 2D ferromagnet. Therefore, 



the exponential temperature dependence of the suscep- 
tibility is the universal property. Certainly one cannot 
seriously rely on the numerical prefactor, the prefactor is 
not universal. Eq. (|10[) gives the magnetic susceptibility 
per lattice site. To compare with experimental data, we 
need to rewrite it per unit area and account for the Bohr 
magneton /is and ^-factor of electrons in silicon, g = 2. 
This gives in units of \ib per tesla per cm 2 



X = ° A1 Nj e 



n lS „7.6JS' 2 /T 



(11) 



where J is taken in Kelvin, ri[ is the density of localised 
electrons per cm 2 and N is the number of electrons per 
droplet. 

The low temperature (T <C JS) and zero magnetic 
field specific heat per lattice site immediately follows 
from the dispersion @ 



C = 



T 



xdx i T , 

To To ■ (12) 



2irJS J e x -l 12 JS 



Hence per unit area the specific heat 



is 



r - —V±2L 

~ 12 N JS 



(13) 



Interestingly, the temperature dependence is the same 
as that for a 2D Fermi gas, Cf = however, the 

behaviour in an external magnetic field is distinctly dif- 
ferent from that of a gas. The specific heat of the fer- 
romagnet is suppressed in the field B ~ T, while the 
specific heat of the Fermi gas is not very sensitive to the 
field. 



III. COMPARISON WITH EXPERIMENTAL 
MAGNETIC SUSCEPTIBILITY DATA 

Experimental susceptibility datai are presented in 
Fig. |21 The MIT point occurs at n = n M iT = 0.85 
(Hereafter for densities we use the units 10 11 cm -2 . Fit 
of experimental data by Eq. (fTTj) gives JS 2 ~ 0.6 K 
for n = 0.55, 0.88 and 1.43, while for n = 4.3 it gives 
JS 2 ps 0.4 K. The fitted curves are shown in Fig. [2] by 
solid lines. For T < 3 K, the fits are very good, while at 
higher temperatures the fitted lines are above experimen- 
tal points. There are two reasons for this, (i) We have 
pointed out above that at higher temperatures Eq. (fTTj) 
overestimates the susceptibility, (ii) It is likely that some 
electron droplets are thermally depopulated with rising 
temperature, reducing the susceptibility. 

The above analysis, in our opinion, is very reliable and 
simple because it is based on a very steep exponential 
dependence. After fixing the exponent JS 2 , we discuss 
now the prefactor in Eq. (fTTj) and here we need to at- 
tract additional arguments. At the MIT and lower den- 
sities, all electrons are localised, hence ni = n. There- 
fore, according to Eq. (HE]), x/rc(T = 1.7 K) - 1.6S 3 /N. 
This should be compared with the experimental value 
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T(K) 

FIG. 2: Experimental datap- for xl n versus temperature for 
various values of electron density n. Below values of n are 
given in units 10 11 /cm 2 . Black circles correspond to n = 
0.55, red squares correspond to n = 0.88, green diamonds 
correspond to n = 1.43 and blue triangles correspond to n = 
4.3. Lines show the exponential fits described in the text. 



x/n(T = 1.7 K) « 3 presented in Fig. [5] Very impor- 
tantly, the prefactor is of the right order of magnitude 
if S ~ N ~ 1. For S = 1/2 and N = 1, the prefactor 
is somewhat small compared to experiment, but still ac- 
ceptable having in mind that theoretically the prefactor 
in Eq. pip is not very precise, so a disagreement by fac- 
tor of several is quite possible. Other combinations with 
S ~ N ~ 1 are also acceptable from this point of view. 
In our opinion, the combination S = 1 and N = 4 is the 
most likely. There are the following arguments in favor of 
this combination, (i) The magnetic susceptibility prefac- 
tor is close to the experimental one. (ii) The experimental 
data o£ii22 show that the saturation spin magnetisation 
in the insulating phase is roughly ^50% of the maximum 
possible spin magnetisation. For S — 1 and N = 4, one 
expects 50% saturation magnetisation, which is close to 
the observed one. (iii) Many-body calculations for small 
2D quantum dots^ show that the ground state spin of a 
four electron dot is S — 1, already at a rather moderate 
value of Coulomb repulsion. Note that in referring point 
(iii), we imply that the two valley dispersion degeneracy 
in 2D silicon is not relevant, the valley hybridisation is 
sufficiently large, at least 3-4 K. Due to the hybridisa- 
tion, electrons occupy only the bonding combination of 
valleys and the droplet is similar to the N = 4 quantum 
dot in GaAs^. Without such hybridisation, the spin of 
the N = 4 quantum dot in silicon would be S = 0, be- 
cause two electrons with total spin S = would occupy 
one valley and the other two electrons also with 5 = 
would occupy the other valley. The n = 4.3 curve in 
Fig. [2] lies significantly below the others. It is natural 
to assume that in the metallic phase (n > Umit)\ the 
number of localised electrons is roughly independent of 
n and is equal to n; ~ timit- With account of this argu- 
ment, the n = 4.3 curve has to be scaled up by a factor 
of 4.3/0.85 « 5. After this scaling, the curve is very close 
to the others. 



All in all, the experimental data are remarkably consis- 
tent with the Heisenberg model picture of electron spin 
droplets. 



IV. HOW TO FURTHER CHECK THE 
HEISENBERG MODEL PICTURE? 

The low temperature exponential divergence of the 
magnetic susceptibility in the 2D ferromagnetic Heisen- 
berg model is necessarily accompanied by high sensitivity 
to the magnetic field. The dependence is clearly demon- 
strated in Fig. [T] where it manifests itself as T m (B) with 
two distinct behaviours for T < T m and T > T m . Al- 
ternatively, one can consider the susceptibility (or mag- 
netisation) as a function of the magnetic field at a fixed 
temperature. This will also have two distinct regimes, 
one with approximately linear dependence of the mag- 
netisation on the magnetic field when B < B* and one 
with very slow increase/saturation of magnetisation at 
B > B* . There are indications of such behaviour in ex- 
isting experimental datai£. Further measurements and 
comparison with results of Monte Carlo simulations can 
shed more light on this problem. 

Another possibility to test the model is to measure the 
specific heat in a magnetic field in the insulating phase. 
We already pointed out in Section II that the magnetic 
field significantly and predictably modifies the specific 
heat. Again, a comparison with the results of Monte 
Carlo simulations for the Heisenberg model would be very 
useful. 

It is known that there are two mechanisms for interac- 
tion between localised spins: (i) usual exchange and (ii) 
superexchange. The exchange mechanism leads to the 
ferromagnetic interaction, while superexchange leads to 
the antiferromagnetic one. Let us denote by R the ra- 
dius of localisation (radius of the droplet) and by I the 
separation between droplets. Obviously R < I. Usu- 
ally at R <~ I, exchange wins and this, in our opinion, 
describes the present situation. On the other hand, at 
R <C I, superexchange always wins^, leading to the an- 
tiferromagnetic interaction. This implies that deeply in 
the insulating phase when ]?<Lwe expect a transition 
to the antiferromagnetic interaction. Since the droplet 
positions are random, the antiferromagnetic interaction 
implies a spin glass state. So, we predict a transition to 
a spin glass state at a sufficiently low density and this 
must be clearly seen in the magnetic susceptibility. Un- 
fortunately we cannot quantitatively predict the critical 
density for the onset of the spin glass state. 



V. CONCLUSIONS 

We suggest the quantum Heisenberg ferromagnct 
model to explain the anomalous magnetic properties ob- 
served in the vicinity of the metal insulator transition in 
a low-density two-dimensional silicon-based electron gas. 
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The ferromagnet is composed of electron spin droplets. 
The observed very steep temperature dependence of the 
magnetic susceptibility is associated with the exponential 
divergence in the Heisenberg model. By comparing the 
experimental data with known analytical and numerical 
results for the 2D Heisenberg ferromagnet, we determine 
the parameters of the model JS 2 « 0.6 K, where 5 is the 
spin of the droplet and J is the ferromagnetic exchange 
constant between droplets. We further argue that most 
likely 5 = 1 with four electrons occupying each droplet 
on average. The 2D Heisenberg ferromagnet is strongly 



and distinctly influenced by magnetic field. Based on 
these properties, we suggest further experiments to test 
the model. 
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